Gyroscopic motion of superfluid trapped atomic condensates 



Halvor M0ll Nilsen f, Dermot McPeake J* and J F McCann f 

f Institute of Physics, University of Bergen, Allegaten 55, N-5007 Bergen, Norway 

I Dept of Applied Mathematics and Theoretical Physics, Queen's University Belfast, Belfast, 

BT7 INN, Northern Ireland 

* NMRC, University College Cork, Lee Makings, Prospect Row, Cork, Ireland. 
Abstract. 

The gyroscopic motion of a trapped Bose gas containing a vortex is studied. We model the 
system as a classical top, as a superposition of coherent hydrodynamic states, by solution of 
the Bogoliubov equations, and by integration of the time-dependent Gross-Pitaevskii equation. 
The frequency spectrum of Bogoliubov excitations, including quantum frequency shifts, is 
calculated and the quantal precession frequency is found to be consistent with experimental 
results, though a small discrepancy exists. The superfluid precession is found to be well 
described by the classical and hydrodynamic models. However the frequency shifts and 
helical oscillations associated with vortex bending and twisting require a quantal treatment. 
In gyroscopic precession, the vortex excitation modes m = ±1 are the dominant features 
giving a vortex kink or bend, while the m = +2 is found to be the dominant Kelvin wave 
associated with vortex twisting. 
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1. Introduction 

Superfluidity is one of the most dramatic illustrations of the quantum nature of matter. It is 
exemplified in its purest form in Bose-Einstein condensation H) of cold atomic gases. The 
quantisation of angular momentum associated with the phase of the quantum fluid has been 
established through the nucleation of vortices EJ |3 recently extending to the observation 
of large vortex lattice structures with grain boundaries [4 J . Experiments have also revealed 
collective irrotational oscillatory flows in the superfluid phase [5 , 6| that revert to thermal flow 
modes and frequencies above the critical temperature. 

One of the fundamental properties of a pure quantum fluid is the coherence of the 
phase of the wavefunction [7| where different regions of the condensate are phase-locked 
in their motion. For example, in the flow of an condensate past an obstacle, the requirement 
of upstream and downstream continuity in mass transport and phase gradient explains the 
capacity for superfluid flow ID, and the transition to dissipation when local phase slip 
arises. For homogeneous and inhomogeneous condensates, there is conclusive evidence, 
both theoretical and experimental, that the nucleation of vortices is the primary mechanism 
responsible for the onset of drag and dissipation in condensate flow [ 8[|9l[T0||. Phase coherence 
is responsible for the quantisation of circulation (angular momentum) and the capacity of 
the fluid to conserve its state of rotation unless acted upon by a critical external moment. 
Consequently, it was predicted ifTTl IT2l that a Bose-condensed gas could be made to exhibit 
gyroscopic motion. Very recently this phenomenon has been observed experimentally lfT3l 
and measured for the first time. In analogy with a spinning top, slightly displaced from 
equilibrium by an impulsive moment, a slight tilt of the axis of rotation (figure [TJ produces 
precessional and nutational motion. In this paper we show that internal flows past the 
vortex create a Magnus force which results in precession, and that this combined with the 
inhomogeneity in the condensate gives rise to vortex bending and twisting into helical modes 
of oscillation. 

The quantum gyroscope is also a unique diagnostic tool able to probe condensate flow 
velocities, acoustic waves, and pressures. The oscillating and twisting motion of a vortex, as 
we will show, is highly sensitive to density fluctuations and inhomogeneities. In this paper we 
study the dynamics of a single-quantum vortex lying near the axis of a cylindrically symmetric 
trap. In our study we have employed four different methods of analysing the motion: (a) 
the Bogoliubov-de Gennes equations (b) the Gross-Pitaevskii equation, (c) hydrodynamic 
equations, and (d) classical equations. While methods (a) and (c) are essentially linearisations 
of (b), they are extremely important in identifying the underlying modes of excitation. We 
study the collective modes of the system and compare with recent hydrodynamic theory. We 
find that while many features of the gyroscopic motion are described by the hydrodynamic 
model, the vortex motion is significantly different. Firstly, the frequencies of precessional 
motion are shifted from the hydrodynamic predictions, secondly the vortex core has a single 
kink, and thirdly the vortex performs a helical Kelvin wave oscillation of mode m = +2 not 
previously seen. 
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Figure 1. Classical rigid gyroscopic motion of the inhomogeneous trapped condensate. 
The figure indicates a surface of constant density for a singly-quantised gyroscope (n = 
1) corresponding to the ground state solution of the Gross-Pitaevskii equation with atom 
interaction strength C = 1000. The classical Euler angles |14|, 9 and <fi, defining the 
orientation of the vortex are shown. 



2. Formulation of the problem 

For a cold weakly-interacting gas, the ground state (condensate mode) dominates the 
collective dynamics of the system. In experimental realizations one can achieve temperatures 
such that T < T c (typically 0.1 to 1/iK) and densities such that the gas is weakly interacting 
and highly dilute. Under such conditions, the condensate of Nq ^> 1 atoms is well 
described by a mean field, or wavefunction, governed by the Gross-Pitaevskii equation, and 
the quasiparticle excitations are acoustic waves within this field. If the perturbations of 
the condensate are small, then it is appropriate and convenient to use the linear response 
approximation, which is equivalent to the Bogoliubov approximation for single-particle 
excitations in highly-condensed quantised Bose gases at zero temperature. 

Consider a dilute system of N atoms, each of mass m a , trapped by an external 
potential V ext (x, t) and interacting weakly through the two-body potential V(x, x'). At low 
temperatures and densities, the atom-atom interaction can be represented perturbatively by 
the s-wave pseudopotential: V(x,x') = {ATxh 2 a s /m a )5^\x — x'), and a s is the s-wave 
scattering length. The dynamics follow from the Hartree variational principle: 



where g = (4irh 2 /m a )N a s , H = —(h 2 /2m a )'V 2 + V ext — fi, and the chemical potential p, 
plays the role of a Lagrange multiplier. The condensate and its excitations can be described 




(1) 
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Figure 2. Surfaces of constant density of the condensate. Results from the solution of the 
Gross-Pitaevskii equation showing vortex bending and twisting, and condensate quadrupole 
oscillations in the scissors mode. A vertical cut through the surface reveals the vortex core 
shape. The images correspond to the parameters C — 1000, A = y7> an d angle 9q = 10°. The 
figure shows a sequence of frames at regular intervals. Top row, left to right t = 0.0, 1.3, 2.6, 
then bottom row, left to right t = 3.9, 5.2 and t = 6.5 in units of u>q . The vertical z-axis 
points in the direction of the axisymmetric trap axis. The unit of length on the figures is 
(n/2m Q ^o) 1/2 . 

by the linear response ansatz: 

4>(x,t) = a (t)<f>(x) + [aj(t)u j (x)e^ t + a*(t>*(cc)e + ^*] (2) 

j>0 

where <\> represents the highly-occupied condensate; that is, \ao\ « vN ^> \aj\, j > 0. From 
the variation 5 (ft*, and linear expansion in the small parameters a J; a* taken as constant, the 
stationary Gross-Pitaevskii equation and Bogoliubov equations follow: 

H <P + g\<P\ 2 <P = (3) 

with ((p, (p) = 1 and J d 3 x f(x)*g(x) = (f,g). The Bogoliubov modes are solutions of the 
coupled linear equations: 

(H + 2g\(f)\ 2 )uj + g(p 2 vj = +h J j j u j (4) 

(H + 2g\<f)\ 2 )v j + g(f)* 2 u 3 = - frw jVj (5) 

Time-reversal symmetry of equations (14151) is reflected in the fact that every set of solutions 
{ujj,Uj,Vj} has a corresponding set {—ujj,v*,u*} and the normalisation can be chosen 
conveniently, such that: (m,Uj) — (vi,Vj) = 5{j. 

2.1. Superfluid dynamics 

The equivalent fluid motion is expressed through the coupled differential equations for mass 
and momentum transport (HI- The mass density p and momentum current density J are 
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Figure 3. Quadmpole oscillations in the vertical plane. Results for (yz)(t) = (2/No)Q yz (t) 
for 9q = 4°, C — 1000 and A = v7 are shown. Gross-Pitaevskii simulation (solid line), 
coherent hydrodynamic state model (dotted line), and the classical model (dashed line). The 
quantum frequency shift leads to a slower precessional motion than predicted by classical 
theory. 
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Figure 4. Spectral analysis P yz (ui) (equation I29> of the vertical -plane quadmpole (yz)(t) 
shown in figure [3] The parameters are 9o — 4°, C = 1000, A = \/7. Results of Gross- 
Pitaevskii simulation (solid line) from equation d28i . coherent hydrodynamic states (dotted 
line) from equation (I39i . Also shown is the classical model data (dashed line) obtained 
from equations ( I29i and d53l >. The higher frequency peak (1) corresponds to the quadmpole 
(h; 0, 1, !)/(«; 0, 2, 1) and the lower frequency (2) to the quadmpole (h; 1, —1, !)/(«; 1, 0, 1). 
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defined as, p = m a ip*ip and J k = {h/2i)(ip*d k ip — ipd k ip*), where the index k denotes the 
kth vector component, and repeated index summation convention is used. The fluid velocity 
is defined by v k = J k / p, or equivalently in terms of the phase, x> °f me wavefunction, 
v k = (h/m a )dkX- Then mass and momentum conservation are expressed by: 

d tP + d J J j = (6) 
d t J k + djT jk + pd k (V ext /m a ) = , (7) 

where the momentum flux density tensor is defined; 

T jk = pvjVk + \5 jk g(p/m a f - (h/2m a ) 2 pd j d k lnp . (8) 
Finally, the mass and momentum equations (6,7) combine to give the acoustic equation: 

d tt p = didjTij + di [p<9i(V cxt /m a )] (9) 

The pressure is defined as [8|, p = \g(p/m a ) 2 — (h/2m a ) 2 pV 2 lnp . The first term represents 
classical pressure, the second term is the quantum pressure and is discarded in the classical 
(hydrodynamic) approximation. The static ( J k = 0) equilibrium density in the hydrodynamic 
limit follows from ©: 

pf\x) = imjg) [p - V cxt (x, 0)] (10) 

For small flow and density fluctuations about equilibrium, the kinetic energy and quantum 
pressure terms are small, so that: T^- w \g{p/m a ) 2 5ij. Using equation © the hydrodynamic 
acoustic modes given by p(x,t) = p^\x) + p^(x)e" lu)t satisfy the linear eigenvalue 
problem: 

-m a u 2 p^ =d t [[p-y ext (r,0)]^pW] (11) 

Finally, the relation of the torque (r, ) by an external moment to the rate of change of angular 
momentum is given by: 

n(t) = d t e ijk XjJ k dVt = - / e ijk Xj [d t T ik + pd k (V/m a )] dtl (12) 
Jn Jn 

where the fluid occupies the volume f2, and is the Levi-Civita symbol. Since = it 

follows that: 



T i(t) = ~ / tijk P Xjd k (V/m a )dtt (13) 



n 



The conservation of circulation, Kelvin's theorem, and its quantisation means that within the 
fluid, only those moments above a critical value can transfer angular momentum associated 
with vortex creation or destruction. Suppose the trap orientation changes with time, then the 
ellipsoidal potential can be written as the quadratic form: V(t) = \dij(t)xiXj where = aji. 
Then the equation for superfluid motion, that is the absence of torque 73 = 0, is simply a set 
of homogeneous linear equations for the quadrupole moments of the condensate: 

a>ki(t)e ijk / p(x,t) XjXi dVt = (14) 
Jn 

These moments are the key physical parameters for the condensate motion. 
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3. Methods 

3.1. Bogoliubov-de Gennes equations 

Consider the atoms confined by an oblate spheroidal (pancake-shaped) trap: V cxt (x,0) = 
\m a uj 2 {r 2 + X 2 z 2 ), where r is the radial coordinate and z the axial coordinate, with A > 1 
the aspect ratio of the trap. For convenience we use scaled dimensionless units for length, 
time and energy namely: (li/2ro I1 uo) 5 ,Wo _1 an( ^ respectively. We define the interaction 
strength by the dimensionless parameter C = 87rN a s (h/2m a tJo)~K so that C — > 
represents the ideal gas and C — » oo, describes the hydrodynamic limit. Experiments studying 
the scissors and gyroscope modes performed by the Oxford group lfT51[T6l[T7lfr? l employ the 
atom Rb 87 for which a s rs 110 a.u.. So a typical value of C corresponding to a trap frequency 
ujq = 2it x 110 Hz and particle number N ~ 5000 would be C ~ 1000, while oo = 27rx62Hz 
with Nq ~ 19000 is equivalent to C ~ 2870. In the ideal gas limit, the energy of the state 
which we label (i; n r , rig, n z ), has the value: 

E(i; n r , n e , n z ) = 2n r + \n e \ + 1 + (n z + |)A (15) 

with n e = 0, ±1, ±2, . . . and the radial and axial quantum numbers are: n r , n z = 0, 1, . . .. 
The corresponding excitation frequencies, with respect to a vortex state (i;0,K, 0), where 
K — 0, ±1, ±2, . . ., are given by: 

u>(i; n r , rig, n z ) = 2n r + \ng\ — \k\ + n z X (16) 

For finite C, the spectrum of excitations must be determined by numerical solution of 
equations (1314ft and ©. Separating variables gives: 

<P(r,z, ¥ >) = 4> K (r,z)e i ^ (17) 

so that the condensate with circulation k is the solution of the equation: 

-{^ + -^ + ir-2-^)^ + W 2 + AV )<^ + C \^A = (18) 

\or z r or oz A r A J 
We use m = n e — n to denote angular momentum with respect to the condensate, that is the 
helicity of the excitations, and it is convenient to use the labelling (h; q r , m, q z ) appropriate 
for the hydrodynamic limit (C — > oo). Thus the quasiparticle amplitudes: 

Mn r , m ,n z (r, z, if) = Un r ,n z ( r , z)^^^ and u^, m ,„,(r, z, if) = Vn r ,n z ( r , z)e i(m - K)Lp (19) 
with corresponding angular frequency, uj nrim , nz , are solutions of the eigenvalue problem 

C(m + K)Un r ,n z {r, z) +C <f) K (r,z) (20) 

C(m - K)v nr:nz (r, z) + C 0* 2 u nr , nz (r, z) = uJ nr , m ,n z v nr ,n z (r, z) (21) 

where 

C(S) = -(& + 7Tr + &-i) k+ ^ + " V)k + 2Clkl * <22) 
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In our calculations, these two-dimensional equations are discretised by Lagrange meshes ifUfl : 
the radial coordinate is defined at M grid points (r 1; r 2 , . . . , r M ) and the axial coordinate at 
iV points (zi, z 2 , ■ ■ ■ , Zn). Therefore: 

M N 

Ur,z) = ^^^(r fc ,^)A fc - V V 1/2 fk(r) 9 i(z) (23) 
fc=i 1=1 

where /, g are Lagrangian interpolating functions such that 

f:(r)f k (r)2irrdr^\ t 5 lk (24) 

9j(z)gii z ) dz m fijSji (25) 



o 



The Lagrange functions for the r-coordinate are chosen to be generalised Laguerre 
polynomials [18|, scaled to encompass the entire condensate, with typically M = 50 mesh 
points. Hermite polynomials are used in the ^-direction so that 

JV-l 

9i(*) = 'E l Xi(zi)Xi(z) (26) 

1=0 

where Xi( z ) — h N 2 w{z)^Hi{z). and Hi(z) are the Hermite polynomials associated with 
weights w(z) = e~ z2 and normalisation factor h N = 2 N n 1 ^ 2 N. A high degree of accuracy 
was found with only iV = 30 points. The resultant eigenvalue problem was solved using 
Newton's method for equation (fT8t and a standard eigenvalue routine for equations (12012 II) . 
Convergence was established by a combination of grid scaling and number of mesh points, so 
that at least six-figure accuracy was assured for all frequencies (see tabled). 



3.2. Time-dependent Gross-Pitaevskii equation 

The time-dependent Gross-Pitaevskii equation, follows from taking arbitrary variation of ip* 
in©: 

[H + g\i>\ 2 -ifi£ t ] V = (27) 

The direct numerical solution of this equation, without linearisation, in combination with 
spectral analysis can also be used to determine the frequency spectrum [fl~8l and density 
fluctuations of the collective excitations. This can be done efficiently and accurately using 
spectral methods. The ground state of the system is found by evolving equation d27t in 
imaginary time using the split-step Fast-Fourier transform propagator lfl~8l . A simple arbitrary 
trial function is used as the initial state. At each time step in the evolution the fluid circulation 
is imposed by a applying a phase gradient corresponding to k = +1 until the excited states 
diffuse from the system and the density and chemical potential stabilise. In practice this 
normally takes a few trap periods at most. Having found the condensate in this manner, the 
system is then allowed to evolve in real time under any external perturbations using the same 
numerical method. 

For example a sudden rotation of the trap by a small angle will disturb the steady state 
flow of the condensate and create a variety of small amplitude excitations. Under these 
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Figure 5. Quadrupole variations in the horizontal plane. Results for (xy)(t) = (2/No)Q xy (t) 
for angle do — 4°, C — 1000 and A = v7- The full line is the Gross-Pitaevskii quantal 
simulation, the dashed line is the classical model. In contrast to the results in figure[3] the two 
models strongly disagree. 
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Figure 6. Helical mode spectrum: the power spectrum P xy (u>), given by equation J29l > of 
the horizontal-plane quadrupole shown in figure |5] The Kelvin mode, (c) rn = +2, is most 
strongly excited and leads to helical oscillations of the vortex core. The mode (d) rn = —2 is 
more weakly coupled, and the lines (a,b) associated with m = ±1 approximately cancel. 
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conditions, the linearisation © is valid and the mode frequencies and densities correspond 
those of the Bogoliubov-de Gennes equations (I4I5I) . It follows that the expectation value of 
any operator will exhibit beats between the various modes. Since a simple rotation does not 
induce centre-of-mass (dipole) motion, the lowest-order distortion will arise in the quadrupole 
moments: 

Q?.(*) = j ij)\x,t) XiXj i(){x,t) dx (28) 

where ip(x,t) is the solution of equation (|27b . The condensate quadrupole moments in 
the vertical and horizontal planes will give information on the vortex bending and twisting, 
respectively. We present results for the scaled moments: (yz) = (2/N )Q yz (t), in figure ©, 
and (xy) = (2/N )Q xy (t), in figure ©, and the corresponding power spectral densities: 

Pij(u)= / e- wt Qij(t)dt (29) 
Jo 

in figures © and ©, respectively. 
3.3. Coherent hydrodynamic states 

The hydrostatic spectrum, equation (fTTT) . also follows from a simplified version of the the 
Bogoliubov-de Gennes equations [ 19 1 when quantum pressure has been neglected: 

q 2 [-V 2 + q-'iV'q)] h+ - \u)h] = (30) 

[_V 2 + q-^q)] (q 2 hj) - $ujhj = (31) 

where hj = Uj ± Vj, and q 2 is the vortex-free hydrostatic density given by 

2 _ {[p - \ {x 2 + y 2 + A 2 ^ 2 )] if ^ - | (x 2 + y 2 + A 2 z 2 ) > 

g Q if/i _i (x 2 + 2/ 2 + A 2 z 2 )<0 

These equations can be solved in closed analytic form using series expansions llT9l to yield 
expressions for the quasiparticle amplitudes. The scissors mode quadrupole (k — 0, m — ±1) 
for A > 1, with p^ h \x) oc xz, first proposed by Guery-Odelin and Stringari [20| and 
experimentally observed by Marago et al. lfl~5l has the frequency: 

Vl s = Vl + A 2 . (33) 

Quadrupoles corresponding to k = 0, m = ±2 are degenerate with frequency: VL xy = \/2. 
The scissors mode, selected by sudden rotation 6 , introduces a perturbation: 

V w m a u%{l - \ 2 )xz (34) 

In the absence of a vortex (k = 0), coupling to the degenerate m = ±1 states is equal and 
opposite, and no helicity arises. If a vortex is present (k = 1), the coupling is asymmetric 
since the state m = +\m\ rotates with the condensate flow at frequency u + , whereas the mode 
m = —\m\ is counter-rotating at frequency uj~ , see figure [T] The frequency difference was 
calculated analytically using a hydrodynamic model |22| with the result that 
u+ u - _ 2(L Z ) _ (L z ) _ (L z ) 
ynaiy H~ 2z i ^xx ^yy 
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Table 1. Table of mode excitation frequencies for A = \/7 as a function of the interaction 
strength: C = &nNna s (h/2m a LUo)~? . The full spectrum is shown in figure^! The selected 
modes in the table are are: (1) (7i;0,l,l); (3) (h; 0,-1,1); (c) (h; 0,2,0); and (d) (h; 0,-2,0). The 
first entry in each column is the result from the Bogoliubov-de Gennes equations (3-5). The 
numbers in brackets correspond to the time-dependent linear response method of calculation. 
The Bogoliubov-de Gennes results for the frequency of precession are denoted by uj p . All 
angular frequencies are in units of ujq- 
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(1.780) 


0.984 


(1.016) 
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3.169 
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(1.370) 


1.708 


(1.744) 


0.111 


(1.105) 


1000 


0.137 


3.085 


(3.079) 


1.201 


(1.194) 


1.640 


(1.644) 


1.197 


(1.195) 



where I xx is the principal inertia moment about a horizontal axis. The result also agrees with 
the result of Svidzinsky and Fetter [11] which treats the vortex core as a perturbation of the 
hydrostatic Thomas-Fermi density. The same result can be found and understood using a 
classical model of the motion (section l3T4b . 

To obtain the amplitude of precession and nutation, the populations of the excited modes 
are required. Since the tilting is small and sudden, we can use diabatic perturbation theory. A 
small sudden rotation by an angle 9 about the y-axis, with generator L y = (2i)~ 1 (L + — L_), 
gives rise to the function: 

Mt = 0) ~ exp [-itr%L y ] <\) h (36) 

The amplitudes are found by projection to the modes in the rotated frame of reference. From 
equation ©: 

j 

and the coefficients a, (t) will be time-independent. Projection by the hydrodynamic state 
hj = Uj — vj gives the time-independent amplitudes aj 

(hj,h+) a, « {hj,exp [-\h- l 6 (L + - L_)] <f> h ) (37) 

The rotation mixes the initial state with m = ±1, ±2, to order 6 , and 9q, respectively. We 
have overlooked the density distortion of the vortex, though included the angular momentum. 
Since the vortex core distorts the density over a comparatively small volume of the condensate, 
then to a first approximation the density is the axisymmetric Thomas-Fermi distribution, 
equation (l32l . Therefore, as a simple first-order approximation we take 

<p h = C^ q(r,z) e iK * 

where q 2 is given by equation d32l . The time-independent equations (30,31) can be used to 
calculate the quasiparticle functions, and the transition amplitudes, aj given by 071) . can be 
extracted analytically as shown in Ell . In this approximation, the excitation spectrum is the 
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same as the vortex free case, and does not include the degeneracy splitting arising in higher- 
order perturbation theory ll22lfTTl . The quadrupole excitations in the vertical and horizontal 
plane are defined as: 

QyM = J 4>* h (x, t) yz M*, t) dx (38) 

Ql v {t) = J V h {x, t) xy Vft(aj, t) dx (39) 

and can be calculated analytically. The results for this coherent hydrodynamic model are 
presented in figures |3] and |4j 

3.4. Classical model 

We use the convention of Goldstein ff4l for the Euler angles (4>,9,ip), describing the 
orientation of a set of rotating axes (x',y',z') with respect to a space fixed frame (figure 
[l]). The principal moments of inertia of the vortex-free density, equation (132b . are: 



l5\N g 



(40) 



(\ 2 _|_ i) 

Ixx = — 7^2 — Is and Izz = l Is where I s = ^m a N 
In general the Lagrangian density, C a , is given by ll23ll : 

h 2 

C a {p, <f> v , p, V0„, Vp; x- 1) = <f) v p - \p(y<p v ) 2 -pU- ^^(Vp) 2 (41) 

orriaP 

where the potential per unit mass is U = (V cxt /m a ) + (gp/2m 2 ) and v = V<p v . Discarding the 
last term gives the classical interacting fluid [24| with the corresponding Hamiltonian density: 

n c = \pv 2 + + > V + y 2 + X 2 z 2 ) (42) 

At equilibrium the condensate has a hydrostatic (vortex-free) density and angular momentum 
N hn. A small change of orientation of the trap (59) creates density and velocity changes, 
p,v: 

d t v + (g/m a )Vp = , d t p + V-( Po Vv) = (43) 

In the body-fixed frame, and for quadrupole (xz) perturbations only, this leads to an effective 
potential 

5V(9) = ±u 2 (l + X 2 ) Po (x 2 + z 2 )(59) 2 (44) 
Then the Lagrangian in the space-fixed frame has the form, for small 9: 
L(9, <p, 9, j>, t) = \l xx sin 2 9 ^ + § I xx 9 2 + \ I zz {i> + <p cos 9) 2 - \l x Sl\9 2 (45) 
From this Lagrangian, Hamilton's equations are 

( j >= (P*-P?™ ) 9 = ^ (46) 

Ixx sin 9 Ixx 
i Pf (p v - Pip cos 9) cos 9 . 

* = IZ 4^ = ° <47) 

Pe= (P.-P^osO)(p a> S e-p,) _ 
1~~ sin 9 
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In anticipation of the numerical results, consider possible analytic solutions of the equations 
for very small angular displacements 9. The top is aligned vertically (9 = 0) and spinning 
with angular momentum J , thus p^ = J and pe = 0. The possible singularity at 9 = 
in the equation for up can be avoided if p v — p^ cos 9 — > as 9 — > 0. A steady azimuthal 
precessional, (p — 0, independent of 9(t), arises from the case p^ = p v = Jq in which case, 



for cos 9 « 1 — |# 2 , we have: 



. Jo 

If) 7H 

y 21, 



XX 



Under the same approximations, and the fact that 9 = pejl xx . 



Therefore the vortex motion is described by 



0«0 o sin v 7 ^ + * 9^^f-t (49) 

which describes a steady azimuthal precession, frequency tUp = J /2I XX , combined with 
sinusoidal nutation at a frequency slightly higher than the vortex-free scissors oscillation. The 
precessional motion splits the frequencies. It is equivalent to a steady background rotation of 
the condensate. In analogy to a Doppler shift, internal flow around the z-axis associated with 
the quantum number ±m would be shifted so that cu + — uj~ = 2u p = Jq/I xx . This splitting, 
understood in terms of classical precessional motion, is identical to the quantum expression 

ifmi22i. 

Initially the condensate has a pure rotation about the z'-axis so that tp = Jo/h, and 
the condensate is impulsively tilted by an angle 9 . The Hamilton equations are solved 
numerically, with the initial conditions, t — 0: 

9 = 9 cp = (50) 

■0 = P v = Jo cos 6*0 
Pe = p^ = Jo 

The time dependence of the products of inertia is governed by the orientation of the axis 
of the top and thus will oscillate in time. The transformations between the laboratory and 
fixed-frame quadrupole gives 

QIM = W*x - hz) sm2# cosy? (51) 

QIM = W** - T ") Sill2d S[Ul P ( 52 ) 

Q c xy (t) = \(I*x ~ hz) sin2y; sin 2 9 (53) 

Since the inertia moments, I xx and I zz , are constant, the time dependence is entirely governed 
by 4>(t) and 9{t). These functions are determined by numerical integration of equations d48t 
with initial conditions (BTT) . The physical meaning of the quadrupole is clear if one visualises 
a rigid vortex motion which nutates and precesses. At the instant the vortex is aligned along 
the y = plane, then (yz) = 0. This happens twice each precession cycle. During the slow 
precession, the vortex line nutates rapidly between the 1st and 2nd quadrants of the yz-plane 




Figure 7. The excitation spectrum of the vortex condensate as a function of interaction 
strength C = %nN^a s (h/2m a LOo)~2 . The solid lines labelled (1), (2) and (3) correspond 
to those found in the power spectrum of (yz). The dashed lines (a),(b) , (c) and (d) are 
prominent in the quadrupole (xy) . The centre-of-mass dipole mode (dotted line) is included 
for completeness as well as the negative frequency mode (h;Q, — 1,0); line (e). The limit 
C — > corresponds to the ideal gas modes (i;n r ,ng,n z ) with frequencies given by equation 
(I16> . Values of selected mode frequencies are given in tabled 



which leads to a change of sign in Q xz each half-cycle of the scissors period. The amplitude 
of the vertical plane (yz) quadrupole is more easily observable by experiment lfT31l since it is 
of order 9 (figure |3]). The horizontal plane motion, on the other hand, is proportional to 6>q 
(figure EJl. 

4. Results and Discussion 

4.1. Excitation modes and spectrum 

Our numerical studies focus on a singly-quantised vortex (k = 1) condensate in a pancake 
shaped trap A = y/1, for detailed theoretical discussion, and A = y/8 ~ 2.83 for comparison 
with experiment. The frequencies of the quadrupole moments can be obtained from the 
excitation spectrum of the normal modes of the system. The spectra obtained for k = 1 
and < C < 1000 are shown in figure |7J and given in table [T] For C ^> 1000 the variation 
in excitation spectrum is slow and smooth, and tends towards asymptotic limits. Consider 
C = 1000 with an equivalent chemical potential for the hydrostatic density in the absence 

2 

of a vortex: /ztf = (15AC/647r)5 « 8.282 3> 1. While this lies well within the region of 
validity of the hydrodynamic approximation, in the sense that ji 3> ujq, the trap energy and 
vortex core inhomogeneity are still significant and one would expect quantal deviations from 
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the hydrodynamic model. Prior to analysing the motion in detail, an example of the density 
fluctuations produced by our simulations is shown in figure |2j The time-dependent Gross- 
Pitaevskii equation d27l) is solved in imaginary time with k = 1 to generate the initial state of 
the condensate. The trap is then suddenly tilted by a small angle, 8 = 10°, about the y axis 
and the evolution of the condensate follows the time-dependent Gross-Pitaevskii equation. 
Surfaces of constant density of the condensate are shown in figure |2]in the tilted space-fixed 
frame of reference. A vertical cut through the surface reveals the vortex core structure over a 
single trap period. The first image t = is the initial configuration, the subsequent images are 
taken at regular intervals, t = 1.3, 2.6, 3.9, 5.2 and t = 6.5, respectively, in units of uJq 1 . The 
scissors oscillation of the condensate corresponds to nutation, that is quadrupole oscillation in 
the vertical plane. This creates flow past the vortex producing a Magnus force which results 
in precession. However when combined with the inhomogeneity in the condensate, the force 
gives rise to vortex bending and twisting into helical modes of oscillation. The modes of 
oscillation can be most clearly identified by spectral resolution of the mass moments. A 
detailed frequency spectrum of the principal modes for gyroscopic excitation is shown in 
figure|7]as a function of interaction strength C. 

4.2. Vortex precession 

Consider firstly the scissors (yz) or (xz) mode of oscillation; the quadrupole moment results 
are presented in figure|3]and the corresponding frequency spectrum in figurelU Data obtained 
using the direct solution of the time-dependent Gross-Pitaevskii equation are compared 
with the coherent hydrodynamic calculation, and the classical top results. The parameters 
C = 1000 and A = y/7 with iVo atoms and 8 = 4° translate to : I s = (6/7)p, TF N , 
I xx = (16/49)/iTFA r and I zz = (A/7)ixtfNq. Considered as a classical rigid top, the 
condensate, viewed along the ^-direction will exhibit quadrupole oscillations due to nutation. 
If the top precesses then this moment vanishes each time the axis nutates along the xz- 
plane. Using equation d53t . the classical moment (yz) oscillates with a carrier frequency 
ujp = (49/32)/i^p « 0.185co> > and carrier amplitude (yz) max = 2(I ZZ — I xx )9 0.283. The 
classical prediction for the beat period of the (yz) quadrupole, is in fairly good agreement with 
the unperturbed coherent hydrodynamic results cu p = |(3 — V^Vo ~ 0.177co> as shown in 
figure |3] Neither of these results compare well with the more accurate lower-frequency result 
from the Gross-Pitaevskii equation: uj p = 0.137a;o- The lower beat (precession) frequency 
(figure |3j contrasts with higher mode frequencies compared with the hydrodynamic limit, 
as illustrated in figure |U In figure [7] it is clear that, for C = 1000, the curves (1) and (2) 
corresponding to these modes have not reached the asymptotic limit C — > oo. 

Very recently measurements of vortex precession 1T3l in a trap were made under the 
following conditions A w 2.83, lu = 2vr x 62 Hz, and N « 19000 ± 4000 atoms. This 
corresponds to an interaction strength C ~ 2870 ± 600 which is within the hydrodynamic 
regime but should exhibit quantal effects. The classical precession frequency is uj p = 
0.120 ± 0.012a> where the error is due to the uncertainty in iV . A calculation using the 
Gross-Pitaevskii equation gives Wp PE ~ 0.104 ± 0.08a;o- As in the case previously discussed 
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(A = \/7, C = 1000) this rate of precession is slightly lower than the classical/hydrodynamic 
prediction. The equivalent frequencies for the vortex precession are /^ PE ~ 6.5 ± 0.5Hz and 
fp ~ 7.4 ± 0.5Hz. Measurements were performed by imaging the vortex line in the horizontal 
plane and trap tilting along both the xz and yz planes and were in good agreement. Taking 
the average of the results gives the experimental estimate / cxp = 7.8 ± 0.6Hz; a slightly faster 
rate of precession than the quantal prediction, and more in line with hydrodynamic theory. 

While the limitations of zero temperature models, as described in this paper, are well 
known, the hydrodynamic approximation is generally less accurate than the Gross-Pitaevskii 
equation. Moreover, important effects neglected by this model are the viscous effects of the 
thermal fluid, as observed lfi~3l . However, these effects would probably reduce ujp FE further 
and thus increase the gap between experiment and theory. It has been suggested [13] that 
two other effects might have a role in explaining the differences. Firstly, the possibility of 
off-centre precession [ 25 ] and secondly, the influence of additional edges vortices lf26l might 
have significant effects. While the results of the observations and the theory outlined above 
are in very good agreement, there is a considerable margin for improving the simulations to 
study the contribution of these effects and of the influence of thermal damping. 

Considering some of the other modes, we note that this scissors excitation of the vortex 
state weakly populates a low-frequency mode at u = 1.201; on figure[4]and labelled line (3), 
and also marked on figure|7] This is the lowest z-dipole mode (i; 0, 0, l)/(h; 0, —1, 1). The 
variation of frequency of the mode is given in table [T] and drawn in figure [7] Checks on the 
numerical accuracy of the numerical methods were also provided by studying the centre-of- 
mass modes whose frequencies remains constant as the particle number or interaction strength 
C varies; a consequence of the Kohn theorem 11271 l28l . figure UJ We note the labelling of 
axial excitation modes (i; 0,2,1) — > (h; 0, +1, 1) and (i; 1,0,1) — ► (h; 1, — 1, 1), reflects 
the quadrupole symmetry in the hydrodynamic regime. Values of the (h; 0, +1, 1) mode 
frequency are given in table[lj in which the highly accurate Bogoliubov de-Gennes results are 
compared with results using the time-dependent linear response method ifTBl . This confirms 
the accuracy both of the direct time-dependent method and the linear response approximation. 

4.3. Kelvin modes 

Kelvin modes are self-induced helical waves of a vortex line and have a well-established 
dispersion relation for infinite homogeneous fluids [29]. However, in a trapped condensate, 
the finite boundaries means these modes are quantised, and the inhomogeneity of the 
condensate means that the fluid pressure varies along the vortex line. The signature of the 
Kelvin modes is the distortion of the vortex line and can be detected in the gyroscope motion. 
The lowest-order and most strongly-coupled modes, m = ±1, are excited by the small angle 
rotation of the trap. Since the modes m = +1 and m = — 1 are populated in approximately 
equal weights this leads to a kink in the vortex in the vertical plane. The bending and twisting 
of the vortex line is more pronounced at the low-density low-pressure edges of the condensate 
(figure |2j). Analysis of the horizontal plane quadrupole moment (xy) shows periodicity but 
irregularity (figure 13). The two-mode classical theory, equation (l53l) . strongly disagrees with 
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the quantal results shown in figure ©. The amplitude of oscillation is much larger than 
classically predicted, and is associated with the vortex s-shape outward bend. However the 
main feature is the quantal oscillations that indicate twisting, that is bending of the vortex 
projected onto the (xy) plane. The horizontal-plane spectrum (figure |6j) is dominated by two 
pairs of lines labelled (a,b) and (c,d). Figure |7] and tabled indicate the frequency variation, 
with atom number C, of the helical modes. The lines (c) and (d) are the quadrupole pair 
(h; l,+2,0) and (h : 1,-2,0), respectively. In the hydrodynamic limit with k = they 
corresponds to the doubly-degenerate xy-quadrupole with frequency, u = y/2. 

The correlation of the pair (c) and (d) to the ideal gas states (z; 1, 3, 0) and (z; 1, —1, 0), is 
shown in figure|6j The splitting of this quadrupole pair Acu = 0.443co>o is a signature of vortex 
precession, though again it is substantially smaller than the classical prediction Auj p m 0.548 
associated with a rigid vortex crossing the x = and y = planes during precession. The 
dominance of mode (c) over (d) is apparent (figure |6j) and means the vortex has a positive 
helicity of m = +2. Indeed, the strong coupling of m = 2 states with vortex motion has been 
used to measure the angular momentum of the vortex state [26|. The pair of high-frequency 
lines (a) and (b) also contribute strongly to the spectrum. These are the dipole-like modes 
(h\ 1, 1, 0) and (h; 2, —1, 0) , degenerate as C — > with symmetry (z; 1, 2, 0) and (i; 2, 0, 0). 
The pair of lines are of approximately equal weighting and therefore cancel helicities. As 
C — > oo, these lines converge slowly towards the common limit uj = 2.321. Finally the 
spectrum shows two low-intensity lines at uj « 0.274 and uj ~ 1.892. These are not new 
modes but rather vestiges of the (yz) modes; the low frequency line being the beat frequency 
corresponding to the splitting of lines (1) and (2) in figure |H and the higher frequency line, 
the frequency difference of modes (1) and (3). 

5. Conclusions 

In conclusion, we have studied the gyroscopic dynamics of a trapped Bose-Einstein 
condensate containing a vortex. We modelled the system using a classical top description and 
numerically by solution of the time-dependent Gross-Pitaevskii equation. We also compared 
our results with the hydrodynamic approximation. The linear-response equations for the 
system were solved giving the excitation spectra and amplitudes. The superfluid precession 
and nutation of the vortex were found to display quantum frequency shifts. The precession 
frequency was calculated and found to be consistent with recent experiments, though a small 
discrepancy exists. Vortex bending and twisting modes of excitation were observed. We 
found, in the scissors excitation of the vortex state, that vortex bending in the vertical plane 
is associated with the m — ±1 modes of oscillation, while the helical oscillation, vortex 
twisting, is dominated by the m = +2 Kelvin mode. 
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